Experimental Observation of Environment-induced Sudden Death of Entanglement 
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We demonstrate the difference between local, single-particle dynamics and global dynamics of entangled 
quantum systems coupled to independent environments. Using an all-optical experimental setup, we show 
that, while the environment-induced decay of each system is asymptotic, quantum entanglement may suddenly 
disappear. This "sudden death" constitutes yet another distinct and counter-intuitive trait of entanglement. 
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The real- world success of quantum computation HI 0, 0] 
and communication EL IS EL IsL I2L [l^] relies on the longevity 
of entanglement in multi-particle quantum states. In this re- 
spect, the presence of decoherence [11] in communication 
channels and computing devices, which stems from the un- 
avoidable interaction between these systems and the environ- 
ment, presents a considerable obstacle, as it degrades the en- 
tanglement when the particles propagate or the computation 
evolves. Decoherence leads to both local dynamics, associ- 
ated with single-particle dissipation, diffusion, and decay, and 
to global dynamics, which may provoke the eventual disap- 
pearance of entanglement HI El El EH EH]. This phe- 



nomenon, known as "entanglement sudden death" Jl6fl . is 
strikingly different from the single-particle dynamics, which 
occur asymptotically, and as a result has been the focus of 
much recent theoretical work HI El El El El] • We have 
experimentally demonstrated the sudden death of entangle- 
ment of a two-qubit system under the influence of independent 
environments. Our all-optical setup allows for the controlled 
investigation of a variety of dynamical maps that describe fun- 
damental processes in quantum mechanics and quantum infor- 
mation. 

Consider a two-level quantum system S (upper and lower 
states |e) and \g), respectively) under the action of a zero- 
temperature reservoir R. At zero temperature, the reservoir R 
is in the (0)^ (vacuum) state, and the S — R interaction can be 
represented by a quantum map, known as the amplitude decay 
channel 111 |2|] : 
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Under this map, the lower state \g) is not affected, while the 
upper state \e) either decays to \g) with probability p, cre- 
ating one excitation in the environment (state (1)^), or re- 
mains in |e), with probability 1 — p. This would be the situa- 
tion, for instance, in the spontaneous emission of a two-level 
atom. In this case, the state would correspond to one 
photon in the reservoir. Under the Markovian approximation, 
p = 1 — exp (— Ft), that is, the decay probability approaches 
unity exponentially in time. As an initial pure state a\e) + b\g) 
decays, it gets entangled with the environment, gradually los- 
ing its coherence and its purity over time. Complete decay 
only occurs asymptotically in time (p — > 1 when t — > oo), 
when the two-level system is again described by the pure state 



\g). 

Now consider two entangled qubits that decay according 
to the map ([]])• How does the entanglement of the two-qubit 
system evolve? Does it mimic the asymptotic decay of each 
qubit, disappearing at t — > oo, or does it disappear at some 
finite time? This question has been explored theoretically ill 
II [3 01 [12], but up to now there has been no experimental 
investigation of the relation between the global entanglement 
dynamics and the local decay of the constituent subsystems. 

In order to adequately answer these questions, one needs 
a formal definition of entanglement. A convenient measure 
of entanglement for a two-qubit system is the concurrence C, 
introduced by Wootters 11711 . and given by 



C = max{0,A} 



where 



A= \Ai 



(2) 



(3) 



and the quantities are the positive eigenvalues, in decreas- 
ing order, of the matrix 



p(<Ty®<Ty)p* (cr y ®CTy) , 



(4) 



where p is the density matrix of the bipartite system. Here a y 
is the second Pauli matrix and the conjugation occurs in the 
computational basis {|00) , 1 01) , 1 10) , |11)}. C quantifies the 
amount of quantum correlation that is present in the system, 
and can assume values between (only classical correlations) 
and 1 (maximal entanglement). 

For the dynamics given by Eq. (Q]), and an initial state of 
the form |<I>) = \a\ \gg) + \/3\ exp(iS) \ee), the entangle- 
ment decay dynamics depends on the relation between \a\ and 
|/3 1 11511 . Concurrence in this case is given by 



C = max{0,2(l-p)|/3|(H-p|/3|)} . 



(5) 



From this expression, one can see that for \/3\ < \a\, entangle- 
ment disappears only when the individual qubits have com- 
pletely decayed (p = 1), while for \/3\ > \a\, entanglement 
disappears for p = \a//3\ < 1, which corresponds to a fi- 
nite time. This phenomenon has been called "entanglement 
sudden death" W16I1 . Since the concurrence of the initial state 
(p = 0) is C = 2 1 a/3 1, the entanglement dynamics of two 
states with the same initial concurrence can be quite different. 
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Photons are a very useful experimental tool for demonstrat- 
ing these properties and, more generally, for studying quan- 
tum channels like the one given in Eq. (Q]), since the decoher- 
ence mechanisms can be implemented in a controlled manner. 
Let us associate the H and V polarizations of a photon respec- 
tively to the ground and excited states of the two-level system 
S. The reservoir R in turn is represented by two different mo- 
mentum modes of the photon. Fig. [T^ shows a Sagnac-like 
interferometer which implements the amplitude-decay chan- 
nel (Q]) for a single qubit. A photon, initially in the incom- 
ing part of mode a, is split into its horizontal (H) and verti- 
cal (V) polarization components by a polarizing beam split- 
ter (PBS1). Let us ignore the half- wave plates HWP1 and 
HWPC momentarily. The F-polarization component is re- 
flected and propagates through the interferometer in the clock- 
wise direction, and, if unaltered, reflects through PBS1 into 
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FIG. 1: Experimental setup, a Amplitude decay channel for a sin- 
gle photonic qubit. Due to the polarized beam splitter PBS1, the 
H and V components of polarized photons propagate respectively 
along counter-clockwise and clockwise paths within the interferom- 
eter. With half- wave plate HWP1 set at 0°, they are coherently re- 
combined into the outgoing spatial mode a, which represents the 
"vacuum state" of the reservoir. For other angles of HWP1, the 
^/-component undergoes a rotation, corresponding to its probabilis- 
tic "decay" into the ^-component, which PBS1 sends to outgoing 
mode b, representing the "one-excitation state" of the reservoir. Wave 
plates HWP2 and QWP2, together with PBS2, are used for tomog- 
raphy of the polarization state of outgoing modes a and b, which 
are recombined incoherently on PBS2 using HWP3, and sent to the 
same detector. IF is an interference filter, b Amplitude decay chan- 
nel applied to entangled qubits. The entangled state is generated by 
parametric down-conversion in type-I non-linear crystals. 



the outgoing part of mode a. The i^-polarization component 
is transmitted and propagates through the interferometer in the 
counter-clockwise direction and transmits through PBS1, also 
into the outgoing part of mode a. However, the interferom- 
eter is aligned so that the two paths are spatially separated, 
making it possible to manipulate the H and V polarization 
components independently. 

To realize the amplitude decay given in map ©, we use 
HWP1 to rotate the polarization of the V component to 
cos(2<9) \V) + sin(20) \H), where 6 is the angle of HWP1. 
Suppose that an incoming photon is F-polarized. When this 
photon exits the interferometer through PBS 1 , it is transmitted 
into mode b with probability p = sin 2 (20) and reflected into 
mode a with probability cos 2 (26). This evolution can thus 
be described by \V) \a) — > V^P \V) \a) + ^/p \H) \b) . 
Identifying the outgoing modes a and b (which correspond to 
orthogonal spatial modes) as the states of the reservoir with 
zero and one excitation, respectively, this operation is equiva- 
lent to that on the \e) \0) R state in map (Q}. An incoming H- 
polarized photon is left untouched, corresponding to the first 
line in Eq. (Q}. This process is therefore identical to the decay 
of a two-level system. Half- wave plate HWPC, oriented at 0°, 
is used solely to match the lengths of the two optical paths. 
The path lengths are adjusted so that if HWP1 is oriented at 
0°, the polarization state in mode a after the interferometer 
is exactly the same as the input state. Photons in modes a 
and b are then directed to the same quantum state tomography 
(QST) system, composed of a quarter- wave plate QWP2, half- 
wave plate HWP2, and the polarizing beam splitter PBS2, and 
then registered using a single-photon detector, equipped with 
a 10 nm FWHM interference filter, and a 1.5 mm diameter 
aperture. Mode b is recombined incoherently with mode a 
on PBS2, so that both modes can be detected with a single 
detector. This is achieved by assuring that the path length 
difference between modes a and b is greater than the coher- 
ence length of the photons, which is determined by the width 
of the interference filters (~ 0.1mm). The half- wave plate 
HWP3, aligned at 45°, transforms i^-polarized photons into 
V -polarized ones. Since PBS 2 transmits i^-polarization and 
reflects F-polarization, the combination of HWP3 and PBS2 
reflects photons that were originally i^-polarized. This as- 
sures that the QST performed is identical for both modes a 
and b. 

Using the interferometer described above, we studied both 
the decay of a single-qubit and the dynamics of two entangled 
two-level systems interacting with independent amplitude- 
decay reservoirs. The experimental setup is shown in Fig. [I]). 
Polarization-entangled photon pairs with wavelength centered 
around 884 nm were produced using a standard source lfl8ll 
composed of two adjacent type-I LiI03 nonlinear crystals 
pumped by a 441.6 nm c.w. He-Cd laser. One crystal pro- 
duces photon pairs with ^-polarization and the other pro- 
duces pairs with i^-polarization. After propagation and spa- 
tial mode filtering, the H and V modes are spatially indis- 
tinguishable, and a photon pair is described by the pure state 
|$) = \a\ \HH) + \[3\e i5 \VV) with high fidelity. A half- 
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wave plate (HWP) and a quarter- wave plate (QWP), placed in 
the pump beam, allow the control of the coefficients |a| and 
|/3|, and the relative phase S of the state 11811 . 

The decay of a single qubit was investigated experimen- 
tally for both H and ^-polarized photons, by generating states 
\VV) and \HH), and registering coincidence counts, with 
one photon propagating through the interferometer and the 
other serving as a trigger. The coincidence detection win- 
dow (cx 5ns ~ 1.5 m) was larger than the path difference 
between outgoing modes a and b (~ 5 cm). Fig. [2^ shows 
Pv(V), Ph(V), Pv(H) Ph{H) as a function of p, where 
Pj(K) is probability of finding an input if -polarized pho- 
ton in the J state after the interferometer. The linear behav- 
ior in p is characteristic of exponential decay in t, given that 
p = 1 — exp (— Ft). Vertical error bars were obtained by stan- 
dard Monte Carlo simulations 112011 . The horizontal error bars 
represent uncertainty in aligning the waveplates. 

For the investigation of the entanglement dynamics, non- 
maximally entangled states were produced, and each pho- 
ton sent to a separate interferometer, which implemented an 
amplitude-damping reservoir, and then to a QST system. The 
half- wave plates HWP1 and HWP4 were set to the same an- 
gle 0, so that the reservoirs, though independent, acted with 
the same probability p. QST of the output two-photon state 
followed the usual recipe consisting of a set of sixteen co- 
incidence measurements 111911 . We repeated the same proce- 
dure for different values of p, obtaining the tomographic re- 
construction of the output two-photon polarization state in all 
cases. The concurrence was calculated using Eqs. © and ©. 

Figure Efc displays the concurrence, and the quantity A, 
given by Eq. 0, as a function of the decay probability 
p, for two initial states that, although not pure, are very 



close to |<£) 



\HH) 



\e ld \ VV): state I, defined by 



J /3 (triangles), and state II, defined 
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(squares). Tomography of the initial states I and II showed 
them to have the same concurrence (~ 0.8), and similar pu- 
rity (~ 0.91 — 0.97). The theoretical curves were obtained 
by applying map (Q} to the experimentally determined initial 
states, corresponding to p = 0. As before, the vertical error 
bars were obtained by Monte Carlo simulation [20]. For initial 
state I, entanglement disappears asymptotically, and the con- 
currence goes to zero only when both individual systems have 
decayed completely (p = 1). For initial state II, however, the 
entanglement behaves very differently: the concurrence goes 
to zero at a finite time (p < 1). This demonstrates "entangle- 
ment sudden death." 

It is also illustrative to study the purity, defined as tip 2 , 
in function of the decay probability, as shown in Fig. [3] for 
states I and II. We see that in both cases the purity reaches a 
minimum but is restored when p = 1, when all photons have 
"decayed" to the i^-polarization state. 

In order to further illustrate the usefulness of the present 
scheme for studying decoherence of entangled systems, we 
have performed a second experiment studying the action of a 
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FIG. 2: Results for amplitude decay channel, a Experimental ampli- 
tude decay for a single qubit. Pv(V) and Ph(V) are the probabili- 
ties of detecting an input ^-polarized photon in the V and H states, 
respectively. Pv (H) Ph (H) are the the probabilities for an input 
H photon. The points correspond to experimental data, and the lines 
are linear fits, b Entanglement decay as a function of the probability 
p. The squares correspond to experimentally obtained values of A 
for the case \(3\ 2 — 3|a| 2 . The solid line is the theoretical prediction 
of the concurrence for this state, given by Eq. [2] while the dotted 
line shows the value of A, given by Eq. The triangles are experi- 
mental values of A for the case \(3\ 2 = \a\ 2 /3, and the dashed line is 
the theoretical prediction for A and C, which are equivalent for this 
state. 



pair of dephasing reservoirs, described by the map (US] 
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This map represents elastic scattering between atom and reser- 
voir. States |e) and \g) are not changed by the interaction, but 
any coherent superposition of them gets entangled with the 
reservoir. There is no longer decay, but only loss of coherence 
between ground and excited states. The dephasing map can be 
implemented with the same interferometer through the addi- 
tion of an extra HWP at 45° in mode b before the QST system 
(or, equivalently, through the removal of HWP3 and redefi- 
nition of the QST measurements). For the dephasing chan- 
nel, pure states I and II present identical behavior, becoming 
completely disentangled only when p = 1. The concurrence 
(squares) and bipartite purity (triangle) as a function of p for 
the entangled state II is shown in Fig. |U The theoretical pre- 
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decay and dephasing. The experimental setup represents a re- 
liable and simple method for studying decoherence of entan- 
gled systems interacting with controlled reservoirs. 
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FIG. 3: Purity as a function of p for the amplitude damping chan- 
nel. The squares correspond to experimentally obtained values of the 
purity for the case \(3\ 2 — 3|a| 2 , while the solid line is the theoret- 
ical prediction. The triangles are experimental values of the purity 
for the case |/3| 2 = |a| 2 /3, and the dashed line is the corresponding 
theoretical prediction. 



diction for the concurrence disappears slightly before p = 1, 
due to the fact that the initial state is not 100% pure. 

In conclusion, we have reported an experimental demon- 
stration of the sudden disappearance of the entanglement of 
a bipartite system, induced by the interaction with an envi- 
ronment. We have shown that entangled states with the same 
initial concurrence may exhibit, for the same reservoir, either 
an abrupt or an asymptotic disappearance of entanglement, in 
spite of the fact that the constituents of the system always ex- 
hibit an asymptotic decay. We have explicitly demonstrated 
that this behavior also depends on the characteristics of the 
reservoir, through two examples, corresponding to amplitude 
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FIG. 4: Experimental results for the dephasing reservoir. Concur- 
rence (squares) and purity (triangles) are shown for the case \(3\ 2 = 
3\a\ 2 . The solid line is the corresponding theoretical prediction for 
concurrence, given by Eq. (J2}. The dashed line is the theoretical 
prediction for purity, given by trp 2 . The concurrence goes to zero 
asymptotically. 
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